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^ 1 Abstract 



We study the azimuthal asymmetries for quark pair production in proton-nucleus collisions using a 
hybrid approach in which the nucleus is treated in the Color Glass Condensate (CGC) framework 
while the Lipatov approximation is applied on the proton side. Our treatment goes beyond the 
large N c limit. We particularly focus on the so-called correlation limit where the imbalance of 
I the transverse momentum of the quark pair is much smaller than the out-going individual quark 

transverse momenta. In this kinematic region, a matching between the hybrid approach and a 



factorization in terms of transverse momentum dependent parton distributions (TMDs) has been 



found. It is shown which of the various unpolarizcd and linearly polarized gluon TMD distributions 
contribute to cos 20 and cos 40 modulations of quark pair production. 



. 1 Introduction 

in 

. During the past three decades, transverse momentum dependent (TMD) factorization framework has 
been systematically developed (TMDs) [1,2]. While many theoretical tools have been developed and 
were applied to a wide selection of hard scattering processes (see, for example Refs. [3,4] and references 
therein) and much insight was gained, a formalism which would be on the one side complete and in 
all aspects well-defined and on the other side easily applicable to all processes has not yet been found. 
Therefore, one still has often to resort to models (like the one used here) for the practical calculations. 
In fact, it seems as if di-jet production in hadron-hadron collisions [5, 6] is one of the especially 
problematic cases. Recently, however, an effective TMD factorization has been established at small 
x for di-jet production in nucleon- nucleus collisions [7,8]. This was achieved by extrapolating the 
full Color Class Condensate (CGC) result to the so-called correlation limit where the /cj_-imbalance 
of di-jets is much smaller than each of the jet transverse momenta. The key observation was that 
in the correlation limit, two out-going jets stay very close in position space, such that multiple-point 
functions appearing in the CGC formalism collapse into the two point function. The derivative of 
the two point function with respect to the transverse coordinate was related to the various gluon 
TMDs. Following this line of argument one would conclude that the basic building blocks are gluon 
multiple-point functions, namely Wilson line, rather than gluon TMDs. Note that the multiple soft 
gluon interaction between the nucleon and the active partons is neglected in the CGC calculation 
since the background gluon field inside a large nucleus is much stronger than that inside a nucleon. 
This leads to the absence of contributions causing a violation of generalized TMD factorization [6]. 
A comprehensive review covering among others the relation between the CGC formalism and TMD 
factorization can be found in Ref. [9]. 

In this paper, following the same spirit, we study quark pair production in the correlation limit in 
pA collisions, with special focus on the polarized case. Quark pair production in high energy hadron- 
hadron collisions has been investigated in the framework of collinear factorization in [10, 11] and in 



1 



kj_ factorization in [12-14]. Later, a CGC calculation has shown that for quark pair production in 
nucleon-nucleus collisions, the result for k± factorization can be recovered from the CGC one in the 
dilute limit where the gluon densities are not too high, but that this fails in a dense medium [15, 16]. 
In addition to the two point function, three point functions and four point functions also show up in 
the CGC calculation. These multiple-point functions provide a deep access into the saturation physics 
in the kinematical region beyond the correlation limit. 

We will first reproduce the full CGC result using a hybrid approach [17] in which the nucleus is 
treated in the Color Glass Condensate(CGC) framework [18, 19] while the Lipatov approximation [12, 
13,20,21] is applied on the proton side. The second step is to extrapolate this result to the correlation 
limit by using a power expansion of the hard coefficients. The fact that the hard coefficients become 
independent of the gluon transverse momentum enables us to integrate out one or two gluon transverse 
momenta. Correspondingly, the three point functions and four point functions collapse and can be 
expressed as derivatives of the two point function. The latter can be related to the various gluon 
TMDs. Using the different polarization tensor structures, the gluon TMDs are classified into an 
unpolarized gluon distribution and the distribution of linearly polarized gluons. The latter one is 
responsible for cos 2(f) and cos 4(f) azimuthal asymmetries for quark pair production in proton-nucleus 
collisions. 

The linearly polarized gluon distribution [22,23] (h^ 9 in the notation of Ref. [24]) recently has 
attracted quite a lot of attention. It is the only spin dependent gluon TMD for an unpolarized nu- 
cleon/nucleus, and may be considered as the counterpart of the quark Boer-Mulders function h^ q (x, k±) 
[25]. However, in contrast to the latter, h^ 9 is time-reversal even implying that initial/final state in- 
teractions are not needed for its existence [26,27]. Despite this fact, it does receive the contributions 
from the initial/final state interaction, leading to the process dependent gauge links. This distribution 
function is of phenomenological interest, especially for small- x physics at RHIC and LHC because a 
calculation in the saturation model [28] showed that its contributions are (at small-x) as large as those 
proportional to the unpolarized gluon distribution. The saturation model calculation also reveals that 
the linearly polarized gluon distributions with different gauge link structures differ significantly at 
low transverse momentum, though they all recover the normal perturbative tail at high transverse 
momentum. 

A few ways of accessing h\ 9 have been put forwarded, namely through measuring azimuthal cos 2(f) 
asymmetries in processes such as jet or quark pair production in electron-nucleon scattering as well 
as nucleon-nucleon scattering [29,30]. Other promising observables are cos 20 asymmetries in photon 
pair production in hadron collisions [31] and in virtual photon-jet production in nucleon nucleus 
collisions [28]. Such measurements should be feasible at RHIC, LHC, and a potential future Electron 
Ion Collider (EIC) [3,4]. More recently, it has been found that the linearly polarized gluon distribution 
may affect the transverse momentum distribution of Higgs bosons produced from gluon fusion [17,32, 
33], color-neutral particles produced in nucleus-nucleus collisions [34], and heavy quarkonium produced 
in hadronic collisions [35]. The authors of Ref. [32] proposed that the effect of linearly polarized 
gluons on the Higgs transverse momentum distribution can even be used, in principle, to determine 
the parity of the Higgs boson experimentally. Transverse momentum dependent factorization has 
been re-examined by taking into account the perturbative gluon-radiation correction to h\ 9 [33] . The 
complete TMD factorization results for Higgs boson production are consistent with earlier findings 
based on the Collins-Soper-Sterman (CSS) formalism [36] and soft-collinear-effective theory [37]. Also, 
the transverse momentum resummation formalism applied to di-photon production in pp collisions [38] 
is closely related. The recent development [39] has shown that it might be also promising to perform 
the resummation procedure on the light-cone. 

The article is organized as follows. In the next section, we start by reviewing our version of the 
hybrid approach which allows us to describe effects caused by the finite gluon transverse momentum 
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on the proton side. Then we reproduce the known result for the quark pair production amplitude 
in pA collisions using this hybrid approach. The next step is a power expansion in the correlation 
limit. We show that the resulting differential cross section depends only on gluon TMDs rather than 
higher multiple-point functions. We particularly focus on the polarized cross section which contains 
the linearly polarized gluon distributions. Our treatment goes beyond the large N c limit used in 
earlier studies. In section III, we discuss our result in the dilute limit, the forward limit and the 
large N c limit. It is shown that our expressions are consistent with the existing results in these 
different limits. In section IV, we re-derive the cross section in the TMD factorization framework. As 
expected, a matching between the CGC formalism and the TMD factorization approach is found in 
the correlation limit. In the end, we summarize our results. 



2 Quark pair production in the CGC framework 

Let us first consider the general case of quark pair production, 

V{P B ) + A{P A /ver nucleon) -»• q(h) + q(l 2 ) + X . (1) 

We assume that the nucleus is moving with a velocity very close to the speed of light into the positive 
z direction, while the proton is moving in the opposite direction. It is convenient to use light-cone 
coordinates for which P£ = Pj^pt 1 and Pg = P^n^ 1 with p = (1,0,0,0) and n = (0,1,0,0). The 
corresponding partonic subprocess is represented by gA{ki)+g p {k2) — > q(h)+q~(h), where k± = x\P^ + 
kyj, denotes the total momentum carried by multiple gluons from the nucleus, and k% = x 2 P^ + k^ T 
is the momentum of the gluon from the proton. In the following the notations k\± and k 2 ±_ are used 
to denote three dimension vectors with kf ± = —k\ T and k 2 ± = —k 2T . To simplify the calculation, 
we choose to work in the light-cone gauge of the proton (A~ = 0). Correspondingly, the polarization 
tensor of a gluon carrying the momentum / is given by, 

e^(l) = -g^ + E ±1 , (2) 
p • I 

As mentioned in the introduction, to facilitate our calculation, a hybrid approach [17] has been 
adopted, in which the nucleus is treated in the CGC model, while on the side of the dilute projectile 
proton one makes the so-called Lipatov approximation [12,13,20,21]. At small x, the gluon radiation 
cascade shows a strong ordering in rapidity. In other words, the radiating color source carries a 
much larger longitudinal momentum than the radiated gluon. It has been shown that a fast moving 
color source can be treated as eikonal line in the strongly rapidity ordered region. Making such a 
replacement is refered to as Lipatov approximation [13]. 

For the process of gluon production in pA collisions, the relevant eikonal line is the past-pointing 
one which is built up through initial state interactions between the color sources inside the proton and 
the background gluon field. The interaction between the classical gluon field and the final state gluon 
emitted from the color source inside the proton does not change this general statement because the 
imaginary part of the scattering amplitude cancels between the different cut diagrams once the final 
states are integrated out. The prescription to treat the eikonal propagator is fixed by this choice. The 
relevant Feynman rules, illustrated in Fig. 1, were given in Ref. [13]. Note that the prescription for 
past-pointing eikonal propagators differs from that for future-pointing eikonal lines. 

The multiple scattering between the outgoing quark pair and the classical color field of the nucleus 
can be readily resummed to all orders [40,41]. This gives rise to a path-ordered gauge factor along 
the straight line that extends in x~ from minus infinity to plus infinity. More precisely, for a quark 
with incoming momentum I and outgoing momentum I + k, the path-ordered gauge factor reads, 

27ra(fc-)p"[tf-l](fc L ), (3) 
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Figure 1: Feynman rules for the eikonal line, which is represented by a double line, a, b and c denote 
color indices. 



with 

-j2„ „— ik i -x 



[U - l](fci) = J d z x ± e- l ^[U(x ± ) - 1] , (4) 
and 

TTf \ 1T> ~*9s f dx~A+(x~, Xx)\ /r\ 

U{x±) = (Ve J -°° ') A , (5) 

where = A+t c with i c being the generators in the fundamental representation. Similarly, multiple 
scattering between incoming gluon (or eikonal line) and classical color field of the nucleus also can be 
resummed to all orders, 

2it5{k~)p* 1 {lJ - l](k x ), (6) 

with 

[U - l](k ± ) = J d 2 x ± e- ik ^[U(x ± ) - 1] , (7) 

and 

U( x± ) = {V e- l9s ^ dx ~ A+{x ~' X±) ) A , (8) 

where (^4 + ) a t> = A£(—if ) with f abc being the generators in the adjoint representation. 

We use these as the building blocks to compute the amplitude for quark pair production in high en- 
ergy pA collisions. It is straightforward to obtain the production amplitude for diagram(a) illustrated 
in the Fig. 2, 

m« = -ig s u(h)yt a - h ~t l ~ m 2 ^. i>[u ] - mh±Hhf-^^n^M^k 2± ), ( 9 ) 

('2 — k\) ~ m + ie k 2 + le p-n 

where the factor 27r<5(£;j~) is suppressed. k 2 = h + h — k\ denotes the momentum of the gluon from the 
proton with l\ and l 2 being the quark and anti-quark momentum, respectively. 4> p (x 2 , k 2 ±) represents 
the probability amplitude for finding a gluon carrying a certain momentum inside the proton, with 
X2d(x2, k 2 ±) = 4> p (j)*. The other diagrams shown in Fig.2 give, 

m* - »l$>ti-Kkw-m^ - _\"^ + )f_ + m T +fc ^ 

t , ,e pa {k 2 ) a k 2 -p 

M® = -i 9 su(hh P t b v(h h 1 ■ ' PCT i fc , 2 | 2 ^U - (fc2 + /£l>P [[/-l] &a (fcix)^(x 2 ,fc 2± ) (12) 

k\ -n — ie (k 2 + k\) z + it p-n 

= -i 9s u(hh p t b v(h) ( ^ + ( %+ % £ qff e ^ pW ^[U - l] ba (h±)M*2, k 2X ) .(13) 
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Figure 2: The diagrams contributing to quark pair production. The gluon line terminated by a 
denotes a classical field insertion. The contributions from all other diagrams disappear because the 
multiple poles are located in the same half plane. 



Putting all these terms together, we obtain the following expression for the complete production 
amplitude, 



M = M {a) + M {b) + M {c) + M {d) +M {e) 



d% 
(2vr) 3 



h - h + Vi + m 



"7 

(li — k\ + k[) 2 — m 2 + it (l 2 — k[) 2 — m 2 + it 

- ig s u(h)Yt b v(l 2 )U ba (k 1± ) 



h - fi - m 



-J 



U(k 1± - k' 1± )t a W(k' 1± )} v^f-^M-n^M^, k 2± 



£p P '(k 2 + ki) e aa i{k 2 ) 
{k 2 + ki) 2 + it kl + it 



k 2 + 



it p-n 



1 £ P a{k 2 + k\) 



k\ -n — it (k 2 + k\) 2 + it 



n a —M x 2,k 2± ). (14) 
p-n 



To arrive the above formula, we have made use of the Dirac equation of motion obeyed by the free 
spinors and the following identity for the Wilson lines in the different representations: 

U(x ± )t a U\x ± ) = t b U ba {x ± ) . (15) 

After some algebraic, this production amplitude can be further rewritten in a more conventional form, 

k 2 -p 



M 



k\ 



-(j> p (x2, k 2 ±) 



xu(h) \ T g (k 1± )t b U ba (k 1± ) + 



d 2 k' 



i± 



(2vr) s 



T qq (k 1L ,k' 1± )U(k 1± - k' l± )t a uHk' 1L )} v{l 2 ) (16) 



where, 



T q g(k 2 ±, kix,k[ ± ) 
Tg(k 2± ,k 1± ) 



dk' dk' + . t U — Iti + $\ + m 

)P 



2tt 



f, 2 - Jjti - m , . , 

TTv? ^—P ( 17 ) 



(h -ki+ k[) 2 -m 2 +it (l 2 - k[) 2 - m 2 + it 



1 



(ki + k 2 ) 2 



kj± ~ + fc 2±) s 
k 2 -p 



k 2 

ki ■ n 



i + 2k' 



2X 



1 



k 2 -p 



(18) 
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This result is in agreement with the production amplitude obtained in Ref. [16] up to a trivial prefactor. 
Note that T qq -(k 2± = 0, k 1± ,k[ ± ) + T g {k 2± = 0, k 1± ) = 0. 

Squaring the amplitude, we obtain the following expression for the pair production cross section, 



da 



dV.S. 



2d 2 k 1± , 2) d 2 k\ , d 2 k'' 1 
■a ho 



(2vr) 4 (2?r) 



2 S 2 (k 1± + /c 2 ± - q±)x 2 g{x 2 , k 2± ) 



x {Tr [(/x + m)T 3? -(/ 2 - m) 7 °7< 7 ] C(x ± ,y ± ,y' ± ,x' ± ) 



k 2 

K 2± 



+Tr 
+Tr 
+Tr 



(/i + m)T 9 ^ 2 - m) 7 Tj' 7 °] C(x ± , y ± , y' ± , y' ± ) 
+ m)T g (} 2 - m) 7 °Tj ?7 ° C(x_l, rr_L, y^, x' ± ) 
(/i + m)T 9 (/ 2 - m) 7 r/ 7 ] C(x ± ,x ± ,yl,yl)} 



(19) 



where g(x 2 ,k 2 ±) denotes the un-integrated gluon distribution of a proton. In the phase space factor 
dV.S. = d 2 h±d 2 l 2 ±dyidy 2 , the quantities yi,y 2 are the rapidities of the produced quark and anti-quark 
respectively. The quark pair imbalance is defined as q± = h±+l 2 ±- The factor 2/(27r) 3 associated with 

phase space integration is chosen such that for a single gluon target, / -j 2 ^~^t^-(U {ki±)W {ki±)) g i uon = 
xi5(l — x\) at lowest non-trivial order (see, for example, Ref. [42]). To obtain the above result, we 
have defined the normalization factor and the flux factor to be k 2 ± /(2k 2 ■ p(N 2 — 1)) and l/(2k 2 ■ p), 
respectively, rather than k 2 ± k 2 ± /(Axix 2 PA • Pb{N 2 — l) 2 ), 1/(4xix 2 Pa ■ Pb) as used in Ref. [13], since 
the Lipatov approximation is only applied on the proton side. We have omitted the arguments of 
T q q,Tg. And T' q q,T' g denote the same quantities with k' 1± replaced by k" ± . The four point function 
C(x_i_,y±,y' ± ,x' ± ) is defined as, 



C(x ± ,y ± ,y' ± ,x' ± ) = Tr c (u(x ± )t a U\y ± )U(y' ± )t a uHx' ± 



(20) 



Here Tr c is a trace over the color indices. The longitudinal momentum fraction of proton and nucleus 
carried by the incoming gluons are constrained by the kinematics, 



X\ 



\hi\er vx + U 2 _i_|e- 



X2 



\h±\e yi +\h ±\e 



in 



(21) 



where y/s is the center of mass energy. 

With this derived full CGC result, we proceed to the correlation limit where \P±\ = \l\± —h±\/^ 3> 
|gj_|/2. In this kinematical region, we may systemically neglect the terms suppressed by powers of 
l&2±|/|-Pj_|> |&:uJ/|-Pl| and Q S /\P±\ in the four hard coefficients. We first perform a Taylor expansion 
of the hard coefficients in terms of k 2 ±. By dropping all terms suppressed by powers of |/c2±|/|-Pl|; 
one ends up with, 



da 
dV.S. 



N 2 -lJ (2tt) 



d 2 k 



J2y £.^11 

d 2 k 2± }± u — S(ki± + k 2 ± - q±)x 2 g(x 2 , k 2 ±) 



d z x±d z y±d 2 x' ± d 2 y' ± e 



(2tt)< 

-ix±-(kix-k' 1± ) „-iy±-k' 1± e ix' x -{kx_L-k'{ x ) Jy' x -k'{ ± 



{ Tr 
-Tr 



(/i + m)f q q{jt 2 - m) 7 °t^ 7 ° 



C{x±,y±,y' ± ,x' ± ) 



(J,i + m)T qq (J 2 - m)j°T^j C(x±,y±,y' ± ,y' ± ) 
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+Tr 
+Tr 



fa + ?n)f g (jl 2 - m) 7 °f^° 



fa+m)f g fa-m) 7 f^ 7 



fc 2 ±=o 



C(x±,x±,y' ± ,x'_ L ) 
C(x±,x ± ,y'_ L ,y'_ L ) \ , 



(22) 



and with T qq , T g given by, 
T q q(ki±,k[ ± ) = i 
f 9 (k l± ) 



dk dk + /i- ^i + ^+m $ 2 ± fe-ft'i-m 

6 ( k l )Pn u I u\2 ^2 i 7771 3 71 £m2 ^2 I „•> ( 23 ) 



2tt 



1 



(fc! + x 2 p B y 



fa -ki + k[) 2 -m 2 + ie (k 2 ■ p) (h - k[) 2 - m 2 + ie 
-2ki± • k 21 _ 



k 2 -p 



-0 + 2$ 



2± 



(24) 



where k 2 ± = k 2 _i/\k 2 ±\ is a unit vector. 

Now let's move on to discuss the power expansion on the nucleus side. The fact that the integration 
over k' 1± , k'{ ± are dominated by the kinematical region k' 1± ~ k'{ ± ~ Q s — because the typical small x 
gluon transverse momentum is characterized by the saturation momentum — allows us to employ the 
power expansion in the correlation limit Q s ~ (ki± + k 2 ±)/2 <C P±. To facilitate the power expansion, 
we replace T qq (ki±, k' 1± ), T g (ki±) with the following two expressions with the help of Ward identities 
(gauge invariance violation terms in the amplitude are proportional to the gluon off-shellness ~ k 2 ± , 
and thus can be neglected in the correlation limit.), 



T q q(k 1± ,k[ ± ) T g q(k 1± ,k' 1± ) 



dk dk 
2^~ 



-8{k{ 



$i - k" x h~h + Vi + m 



C2± 



h 



777 



(xi - x' x )PX fa ~ h + k[) 2 - m 2 + ie (k 2 • p) fa - k[) 2 - m 2 + ie x[P^ 



T g (k l± ) T g (k 1± ) 



x\Pj[ fa — k\) 2 — m 2 + ie (k 2 ■ p) 

h± fa-h- m Ji 

(k 2 ■ p) fa - ki) 2 - m 2 + ie x\P^ 
1 



fc, I =0 



2 + 



fe lx =0 

7.2 
K 1± 



ki± + C( p2~) ~ T gq - i (k\ 1 _ — k 1± ) + T q ^kij_ (25) 



(ki + x 2 P B ) 2 \ x x x 2 P\P B 

1 ~2k 2 ±A ' 



C2± 



(ki+x 2 P B ) 2 k 2 -p 



p —rp.W, 



(26) 



J fci , =o 



with i denoting the transverse index. By making the above replacement, the differential cross section 
can be rewritten in the form, 



da 



dV.S. 



(N 2 - 1) J (27T) 



d 2 ki±d 2 k 21 _ 



x ( Tr 



+Tr 



8 2 (k 1L + k 2± - q±)x 2 g(x 2 ,k 2L ) J d 2 x ± d 2 x' ± e- lk ^< x ^- x '^ 

d 2 C(x ± ,y ± ,y' ± ,x' ± ) 



fa + m)T q lS2 ~ myffgfP 



-Tr 



fa + m)f, f q>i fa - m) 7 °T 1- t ;7 



k2±M±=0 



kn I ,fei | =0 
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-Tr 



-Tr 



-Tr 



-Tr 



-Tr 



+Tr 



(/ 1 + m)f q lS2-mh°f q ^ jl 

(/ 1 + m)r4 i (/ 2 -m)7 Tj; j 7° 
ft + m)T* S2 ~ m)7°C 7 ° 



(i 1+m )T 9)i (/ 2 - m ) T o 5; 7 ° 



*:2X,fcl±=0 



k2±M±=0 



d 2 C(x ± ,y_ L ,y'_ L ,x' A 
dy\dy 3 ± 

d 2 C(x ± ,y ± ,x' ± ,x' ± 



(/ 1 + m)f 3ii (/ 2 -m)7 T g t '. 7 



k2±M±=0 



dx l ± dx'J_ 

d 2 C{x ± ,y J _,x' ± ,x'_ L ) 
dy l ± dx']_ 

d 2 C{x ± ,x±,y' ± ,x'_ L ) 
dx l ± dx'J_ 

d 2 C(x ± ,x ± ,y' ± ,x' ± ) 
dx l ± dy'J_ 

d 2 C(x±,x±,x' ± ,x' ± ) 



J x±=y± 



%±=y± 



c '±=v'± 



dx\ dx 



',1 



(27) 



where the transverse momenta k' 1± and k'-[j_ have been integrated out. As a result, the four point 
functions collapse into the two point functions. The calculation of the Dirac traces in the above 
formula is rather easy, while the evaluation of the soft part in the McLerran-Venugopalan(MV) model 
is a bit more involved. In general, the tensor structure of the soft part can be decomposed in the 
following way, 



d 2 C(x ± ,y±,y'_ L ,x' 1 



dx± idx' 



S -iF 1 (xi,k 1± ) + (k\ ± P 1± -^S l l 



where / denotes / d 2 x±d 2 x'_ L e~ lkl± '( x± ~ x '-i-\ k\ ± is a unit vector k\ ± = ^ij_/| | , and 5j_ = —g l i + 
[p lq n? +p ) ri L )/p ■ n. The four point function C(x±, y±, y'j_, x' ± ) has been evaluated in the MV model 
in Ref. [16]. With the derived four point function, the coefficients Fi and H{ can be computed in a 
tedious but straightforward way. One finds, 



F 3 



d 2 C(x±,y±,y' ± , 


x' ± ) 


dx\dx^_ 




d 2 C(x ± ,y ± ,y' ± , 
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dx 1 , dx 



13 




and, 



Hi 



d 2 C(x±,y±,y' ± ,x'j 
dx± >i dx' Xj 
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with, 
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dx l , dx 
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J-9 



7V2 



h i 9 Dpi x i: h±) - ( 1 - T72 ) ^gC^i) - Tj^i;W( a; i> 



iV2 



27r 4 iV c a s xi 2/i] L fjp(xi,%i_ 



(32) 



(33) 



(34) 

(35) 

(36) 
(37) 
(38) 

(39) 
(40) 



To arrive at the results given above, we have neglected the logarithmic dependence of the saturation 
momentum on r\. Gp>p,Gy\/w > hfop anc ^ hfww are * ne unpolarized gluon dipole distribution, the 
Weizsacker- Williams (WW) type unpolarized gluon distribution, the dipole type linearly polarized 
gluon distribution, and the WW type linearly polarized gluon distribution, respectively. In the MV 
model, they read [28,43,44], 



xiG DP (xi,ki±) = xih^ppixi, fci x 

1 S ± 



CpS± 2 

2AT /Cl± 



d 2 r ± 
(2vr)2 



r 2 Q 2 



x\Gww{xi, ki±) 
x i h iww( x ^ k ^) 



N 2 

c 

N c 4vr 4 a 
N 2 -l 



1 



d 2 -ife 1± -r x / 1 _ 

r. 



8vr- 



S± / dr 



^2(|^i±lk±) 



1 — e * 



(41) 
(42) 
(43) 
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Here S± is the transverse area of the target nucleus. Q 2 = a s N c /jLA^ l/{ r ±^QCD) 1S ^ ne gl uon 
saturation scale with [ia being a common CGC parameter. J 2 is the second order Bessel function. 
Note that our convention for h-^y^w differs from that in Ref. [28] by a factor 1/2. The WW type 
gluon distributions have a clear physical interpretation as the number density of gluons inside a 
hadron/nucleus, while the dipole type distribution does not. On the other hand, the dipole type 
unpolarized gluon distribution in the adjoint representation enters the single gluon production cross 
section in pA collisions [45]. Besides these widely used gluon TMDs, two novel ones are given by, 

d 2 r ± 



XlG q q(xi,kl±) 



2tt 2 



a. 



(27T) 



- ik ^ - r± Q 2 s e 



N 2 



1 



,qq 



8vr 3 



(44) 



(45) 



Collecting all pieces together, the differential cross section for quark pair production can be written 
in the following general form, 

,2 



da 



a 2 s N c 



dV.S. s 2 {N 2 - 1) 



m 



AQ±) + -^^(qi) cos 20 + C(qi) cos - 
-i 1 



(46) 



where <\> is the azimuthal angle between the transverse momenta q± and P±. The coefficients A(q 2 _), 
B(q±) and C(q 2 _) contain convolutions of various gluon TMDs. Instead of presenting the full results 
for these coefficients, we neglect all higher powers in m 2 /P]_ , 



= J d 2 k 1± d 2 k 2± 5 2 (k 11 _ + k 2± - q±)x 2 g(x 2 ,k 2± ) 



[u 2 + i 2 ) 
4ut 



-xiG DP (xi, feu.) + xi 



4 \ 



N 2 



C l-L, 



(47) 



B(q] 



= J d 2 k 1± d 2 k 2± 5 2 (k 1± + k 2 ± - q±)x 2 g(x 2 ,k 21 _) 
{[2(fci 



i± • q±) 
2 



(*-& ) 

s 2 



N 2 ) 



) hi <q q(xi, k\±) + ^ 2 h i ww 



(xi, h±) 



+ 



2(^2± • q± 

(i - u? 



s 2 



xiG DP (xi,ki±) +xi 



4 

N 2 



G q q{xi, k 1± ) + j-^G W w(xi,k 1± ) 



(48) 



C (Q±) = J d 2 k 1± d 2 k 2± 5 2 (k 1± + k 21 _ - q±)x 2 g(x 2 ,k 2± ) 

(2(1? • h±)(q ■ h±) - ki± ■ h±) - 



(i-uf 



-X\l 



ht 9 DP (xi,k 1± ) +xx ( 1 - —5) hi 9 Jxi, ki±) + TTo^tw^i' fe i 



(49) 



where s = (x\Pa + x 2 Pb) 2 , t = (x 2 Pb — h) 2 and u = (x±Pa — h) 2 are kinematical variables defined 
in the usual way. This is the main result of our paper, which opens the possibility to study saturation 
physics through measuring polarization observables as both unpolarized gluon TMDs and linearly 
polarized gluon TMDs are unambiguously determined in the MV model. In the end, we would like to 
mention that it is also feasible to take into account small x evolution effect [46,47]. 
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3 Dilute limit, forward limit, large N c limit 



In the correlation limit, the low gluon densities limit is reached in the kinematic region Q 2 <C kf ± <C 
Pj_. When Q 2 <C k 2 ± , all six gluon distribution functions become identical, though they differ 
significantly at low k±±, 

xiG(xi, k lx ) = Xih^pixi, k 1± ) = xih r $y W (x\, kxx) = xih^ixi, k 1± ) 



xiG DP (xi,ki±) = xxGww{x\,kii_) = xiG q q(xi,ki±) ~ S±— £-3 



N 2 - 1 ha 



k 2 



(50) 



i± 



Therefore, we have, 
A(ql 



J d 2 ki ± d 2 k 2 ±5 2 (ki± + k 2± - q±)x 2 g(x 2 , k 2 ±)xiG(x 1 , ki±) 



N 2 _ I + t 2 t 2 + U 2 



2N 2 



fit 



s 2 



(51) 



y d 2 /ci ± d 2 /c 2 ±5 2 (fci± + & 2 _L ~ Q±) x 29(x 2 , k 2 ±)x 1 G(x 1 , k i± ) 



x4 



iV 2 



2iV2 



2(fcii. • <7±) 2 + 2(k 2± ■ q ± ) 



(52) 



J d 2 fci^d 2 /c 2 ±^ 2 (fci± + &2± - q±)x 2 g(x 2 ,k 2± )xiG(x 1 ,k 

2(9 • fcl±)(? • &2±) - fel± • 



X4 



iv 2 -: 

2iV2 



£2 



1± 
2 1 



(53) 



Here the known unpolarized Born cross section for qq production through gluon fusion has been 
recovered for the unpolarized term, as it should be. Agreement is also found between our results and the 
explicit expressions of the polarized cross section given in the papers [29,30], provided that these results 
are extended to the small x region and the same dilute limit is taken. Note that one automatically 
takes into account the linearly polarized gluons inside a proton in the Lipatov approximation. By 
observing these consistences, we conclude that in the dilute limit, the contribution to gluon TMDs 
from initial/final state interactions encoded in the gauge links can be neglected, and single gluon 
exchange dominates the processes. 

Let us now discuss the expression in the forward limit. Since k 2 ± is relatively small in the forward 
limit, we may make the approximation 5 2 (ki± + k 2 j_ — q±) « S 2 (ki± — qj_) and integrate out ki± and 
k 2 ±- This leads to the following simplified result, 



A(qj_) = x 2 g(x 2 ) 



(u 2 + i 2 ) 



^ — xiG D p{xi, q±) + x\ 



N 2 



Gqq(xi,q±) + 1 r7yGww(xi,q± 



N 2 



B(q±) = x 2 g{x 2 ) 
.J(*-«) 2 



%ihi 9 DP (x 1 ,q±) + 



x\ 



4 \ 



,qq 



1 ~n 2 ) 9-J-) + -mKww( x i^± 



±9 



N 2 

c 



C(q 2 ± ) = , 



(54) 

(55) 
(56) 
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(d) 




Figure 3: The diagrams contributing to quark pair production in TMD factorization approach. The 
gauge link structure of gluon TMD distributions associated with each diagram are different. The 
mirror diagrams are not shown. 

where g{x2) is the integrated gluon distribution function for the proton. One finds that the cos 40 
modulation drops out in the forward limit. In order to compare with the existing result for quark pair 
production in pA collisions in the forward limit, we further take the large N c limit, 



A{ql 

c(«i) 



, {u 2 + i 2 ) \ {i-uf 

x 2 g{x2) — -Tf — < ^ — xiG D p(xi,q±) + x 1 G q g(xi,q±) 

4ut { s z 

x 2 g(x 2 )S^^^ 9 



x 1 h^ 9 DP (x 1 ,q±) +xih^ 9 qq (x 1 ,q ± 



s 2 



. 



(57) 

(58) 
(59) 



where the unpolarized cross section is in agreement with that obtained in Ref. [8] if one uses the 
relations x\Gdp{xi, q±) = Fg l g — Tg 2 ) and x\G q q(x\, q±) = Fgg + Fgg valid in the large N c limit. At 
this point, we would like to emphasize that the large N c limit is not necessarily a good approximation. 
In particular, N c suppressed terms are the dominant contribution in the dense medium since the 
various gluon TMDs scale as: x\ Gp>p k\±) = x\hf 9 DP (xi,ki±) ~ k 2 ± /Q 2 , xiG qq -(xi,ki±) ~ 1, 

x 1 G ww {xi, ki±) ~ ln(Q 2 /k 2 ± ) and xih^ 9 ww (xi,k li _) ~ xihAJx^k^) ~ fiA/Q 2 at low k x 



cij_. 



4 Quark pair production in TMD factorization 

For quark pair production in hadronic collisions, generalized TMD factorization is not valid any longer 
as shown by presenting a direct counter-example in Ref. [6]. However, in pA collisions, if one only 
takes into account the interaction between the active partons and background gluon field inside a large 
nucleus while neglecting the longitudinal gluons attached to the proton side, the type of graph(for 
example, please see Fig. 11 in the paper [6]) which can produce a violation of generalized TMD- 
factorization disappears. In this section, we use this approximation. Admittedly we cannot quantify 
the systematic errors introduced by it. After neglecting the extra gluon attachment on the proton 
side, the multiple gluon re-scattering between the hard part and the nucleus can be resummed to all 
orders in the form of a process dependent gauge link. Due to the different color structures, the gluon 
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TMDs associated with different Feynman diagrams correspond to different gauge link structures. For 
example, the gluon TMD correlation function associated with graph Fig. 3(a) takes the form [48], 



9, 0) 



(2vr)3F+ 
x(P|IV c ^ J P i (0 



N, 



2 Tr 



c/Ot 



N? - 1 N r 



C/Ht 



1 



N2-1 



C/Wt 



F j (0)uW \ \P) , (60) 



where i, j denote the gluon polarization index. The gauge links J/W, 1 are defined as, 



Ul+] 



Ve 



Ve 



ig B f™d<;-A+{<;-,o ± ) Ve -ig s fil dC-i4+(C-,£j.) 



J" 00 dc- A+(C-,o x ) p -iff. dC-A+(C- 



(61) 
(62) 



And C/Pl = C/Wf/Ht = [/[-lt[/[+] emerges as a Wilson loop. At small x, this gluon TMD can be 
expressed as the derivative of a multiple-point function and subsequently be computed in the MV 
model. In order to derive this relation, we made use of the Fierz identities, 

Cix^y^M = Tr c (U(x ± )t a U\y ± )U(y' L )t a U^x' ± )) 

= ^Tr c (uHx' ± )U(x ± ))Tr c (U^y 1 _)U(y , ± )) - ^Tr c (U(x ± )tf(y ± )U(y' ± )tf(x' ± )), (63) 



and the formula, 



/oo 
dx~U[— oo, x~, x±]d l A + (x~ , x±)U[x~ , oo, x±] 
-oo 



(64) 



With the help of the above two identities, one finds, 



2N C 2 f d 2 x ± d 2 x' 



(») m - l 



f-V.s 



l p -ifcu(ii-ii) 



(2tt)' 



5 2 



dxj_ ti dx' ± j 



—C{x_ L ,y_ L ,y ± ,x_ L ) 



(65) 



The normalization on the right hand side of the equation is fixed according to the arguments made 
in Ref. [8]. Following a similar procedure, for the gluon distribution correlation functions associated 
with other diagrams shown in Fig. 3, we obtain, 



9 (b) 



(c) 



(e) 



*(/) 



2A" r 



2A^ C — 



(2vr) 4 



<) 2 



dy±,idy 



—C(x ± ,y ± ,y' ± ,x'_ 



%±=y±, x' ± =y' ± 



d 2 x±d 



(2vr) 



x ± c -ifci i (a i -a/J 



t) 2 



1 2 /• d 2 x ± d 2 x f ± ifc] , (x , _ g > ) 



(2^ 



dx ±4 dy' ±j 

d 2 



—C(x ± ,y±,y' ± ,x'j 



2 2 f d 2 x±d 2 x' ± ^_ ik , , (x , _„/ ) 



(2tt)< 



dx±,idx' ±j 

d 2 



- — C(xj_, a?j_, Xj_) 



2 2 f d 2 x ± d 2 x' 



N r a s 



{2nY 



± e -ik 1± {x ± -x' ± ) 



dx_L :i dx' ±j 
d 2 



—C(x ± ,x ± ,y ± ,xj 



x '±=y'± 



dxxjdy'^ 



—C(x J _,x_ L ,y ± ,x_ 



(66) 

(67) 

(68) 
(69) 

(70) 
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With these relations, all of the unpolarized and linearly polarized gluon TMDs can be calculated in 
the MV model. On the other hand, it is straightforward to compute the partonic hard cross section 
contributions from each diagram in the Fig. 3. Combining the derived gluon TMDs and hard parts 
and summing the contributions from all diagrams, we obtain the finial result in the TMD factorization 
framework. In order to compare the obtained TMD factorization result with that calculated in the 
CGC formalism, we have to take the same dilute limit on the proton side, which means that the 
unpolarized gluon distribution and the linearly polarized gluon distribution inside a proton become 
identical. After making this assumption, a perfect matching between the CGC formalism and TMD 
factorization is found in the correlation limit. Note that this conclusion is valid beyond the large iV c 
limit. 

5 Summary 

In this paper, we have studied quark pair production in high energy proton-nucleus collisions in 
the central rapidity region and in the correlation limit where the total transverse momentum of the 
quark pair (q±) is much smaller than the transverse momenta of the individual quarks (~ P±)- Our 
main focus lay on the polarized case. We first used a hybrid approach to reproduce the full CGC 
result for quark pair production beyond the correlation limit. Our hybrid approach allowed us to take 
into account finite gluon transverse momentum effect on the proton side in a certain approximation. 
Employing a power expansion in the correlation limit, the multiple-point functions appearing in the 
full CGC result collapse into two-point functions and are thus given by a combination of gluon TMDs. 
All finite N c terms are kept in our calculation. The resulting cross section contains cos 2<p and cos 4<p 
dependent terms, where (f> is azimuthal angle between the transverse momenta q± and Pj_. In addition 
to WW and dipole type linearly polarized gluon distributions, the novel linearly polarized gluon 
distribution h±qg also accounts for cos 2^ and cos 4^ modulations. Such observables in principal can 
be measured at RHIC and LHC. 

We further discussed our results in the dilute limit, the forward limit and the large N c limit, 
and found consistency with existing results in the different limits. The physical implications of the 
observed consistences were also addressed. In the end, we showed that a calculation based on TMD 
factorization leads to the same result as obtained in the hybrid approach. The technique introduced 
in this paper can be extended to study di-jet production in the other channels for pA collisions (e.g. 
di-jets initiated by different partons and/or various polarization channels). For these the linearly 
polarized gluon TMDs with different gauge link structures may also manifest themselves through 
cos2<^> or cos4(p azimuthal dependences of the cross sections. We would expect that exploring these 
polarization obervables at small x could open a new path in spin physics as well as in saturation physics. 

Acknowledgments: One of us (Jian Zhou) thanks Andreas Metz for suggesting this work to him 
and for helpful discussions. This work has been supported by BMBF (OR 06RY9191). 
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